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Thickness-mode  quartz  resonators  have  long  been  used  for  timing  and 
frequency  control.  It  is  well  known  that  a  change  in  temperature  results 
in  a  change  in  the  resonance  frequencies  of  such  quartz  plates.  Thus,  a 
major  requirement  for  precision  frequency  applications  is  that  these 
temperature  deviations  are  minimized.  This  report  summarizes  a  detailed 
numerical  investigation  of  the  temperature  sensitivity  of  doubly  rotated 
quartz  resonators.  Various  measures  of  sensitivity  are  used  for  the 
temperature  range  -40  to  8*°C.  Several  cuts  of  quartz  with  stable  tempera¬ 
ture  behavior  are  described  and  compared. 
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Thickness-mode  quartz  resonators  have  long  been  used  for  timing  and 
frequency  control.  It  is  well  known  that  a  change  in  temperature  results 
in  a  change  in  the  resonance  frequencies  of  such  quartz  plates.  Thus,  a 
major  requirement  for  precision  frequency  applications  is  that  these 
temperature  deviations  are  minimized.  This  report  summarizes  a  detailed 
numerical  investigation  of  the  temperature  sensitivity  of  doubly  rotated 
quartz  resonators.  Various  measures  of  sensitivity  are  used  for  the 
temperature  range  -40  to  80°C.  Several  cuts  of  quartz  with  stable  tempera¬ 
ture  behavior  are  described  and  compared. 
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We  begin  by  sumnarizing  the  theory  of  piezoelectric  plate  vibrations. 
We  assume  the  geometry  of  the  plate  depicted  in  Figure  1.  The  plate  is  of 
thickness  2h  and  has  the  y-axis  as  a  normal.  Furthermore,  it  is  assumed  to 
be  infinite  in  extent  in  the  x-z  plane.  These  geometrical  considerations 
imply  that  and  vanish  for  the  field  quantities  and  we  thus  have 
one-dimensional  motion  in  the  y-direction,  i.e.,  the  so-called  simple 
thickness  modes  of  vibration. 


Figure  Is  PLATE  GEOMETRY 


In  the  following,  a  subscripted  index  following  a  coruna  denotes 
differentiation  with  respect  to  that  spatial  index.  To  further  simplify 
the  notation,  we  will  adopt  the  suimation  convention  that  repeated  indices 
are  to  be  sunned  over.  The  stress  equations  of  motion  are: 


m  ++ 

Tij,i  "  p»i 


(1) 


where  T^  is  the  stress  tensor,  p  the  material  density  (assumed  constant) 
and  Mj  is  the  jth  component  of  the  displacement  vector.  Under  the  assump¬ 
tion  of  an  infinitesinal  deformation,  we  have  the  Cauchy  strain  tensor: 

1 

-  (  **i,j  +  (2) 


Sii  " 


The  electromagnetic  characterization  of  the  problem  ;  achieved  using  a 
quasi-static  approximation  which  produces  the  equati  for  the  electric 
displacement  D^,  the  electric  field  Ep  and  potential 


Di,i  =  0 


Ej  *,j 


(3) 

(4) 


The  retaining  equations  are  obtained  from  thermodynamic  considerations. 
Let  the  naterial  under  consideration  be  characterized  by  the  elastic 
stiffness  tensor  at  constant  electric  field,  C^p  the  piezoelectric 
stress  tensor,  e  —  k,  and  the  dielectric  tensor  at  constant  strain  «  ^j. 
These  material  constants  are  related  to  the  field  quantities  by  the  linear 
piezoelectric  constitutive  relations: 


Tij  “  cijklskl  “  ^ij^  ^ 

Di  =  eiklskl  +  ^in^Vi  (6) 

Equations  1-6  can  be  combined,  yielding  the  simplified  4  equations: 

Cijkluk,li  +  ekij  *,ki  =  p  ?j  (7) 

eikluk,li  “  €ik  *,ki  =  0  (8) 

Our  geometrical  considerations  imply  that  <p  p  uk  p  vanish  unless 
l=i=2,  so  7  and  8  become  for  the  infinite  plate: 


C2jk2uk,22  +  e22j  *,22  =  p  "j  (9) 

e2k2uk,22  "  ^22  *,22  *  0  (10> 
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We  now  use  (10)  in  (9),  yielding 


2jk2  + 


*22je2k2 


22 


uk,  22 


=  - p w^u^ 


(11) 


where  u  is  assumed  to  have  a  harmonic  time  dependence.  We  define  the 
piezoelectrically  stiffened  elastic  constants  by: 


c 


2jk2  =  C2jk2  + 
Wien  (11)  becomes  a  wave  equation: 

r 


;22j  e2k2 
€  22 


(12) 


2jk2uk,22  +  pwuj 


+  p  w^u^  =  0 


(13) 


Note  that  this  is  really  3  coupled  equations.  The  solutions  we  seek 
are  acoustic  plane  waves  of  the  form: 


Uj(x,y,z,t)  =  9je 


i  (tt*--wt) 
v 


(14) 


where  ?■  is  a  constant.  Substitution  of  (14)  into  (13)'  yields  the 
equation 


P 


2jk2  ?k  -  pv  3j  *  0 


If  we  define  the  "matrix 
written  as: 


ix"T 


such  that 


r  -  r 

Ijk  '  '2jk2' 


;15) 

then  (15)  mav  be 


r 


P  -  *  9  =  0  (16) 

where  x  =  p  v2  and  ?  is  the  vector  with  components  32,  30(3  ^3*  In 

this  form,  we  clearly  have  an  eigenvalue  problem  (although,  strictly 
speaking,  P  is  not  a  matrix,  i.e.,  a  second  rank  tensor).  The  eigenvalues, 
X(nj,  determine  the  propagation  velocities  of  the  three  modes  by: 

1/2 


V(i)  " 


(-H 


(17) 
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The  solution  of  the  eigenvalue  problem  yields  3  roots,  or  modes,  which  are 
ordered  by  the  relation: 

*  b  *c  (18) 

and  the  convention  is  to  refer  to  the  modes  by  the  letters  a,  b  and  c.  The 
eigenvector  corresponding  to  a  mode  gives  the  direction  oosines  of  the 
displacement. 


In  addition  to  the  eigenvalues  and  eigenvectors,  there  is  another 
parameter  which  is  important  in  characterizing  the  modes.  This  is  the 
piezoelectric  coupling  constant,  which  is  a  measure  of  the  ratio  of 
electrical  energy  to  mechanical  energy  of  the  vibrations.  Following 
Ballato  (1977)  and  Tiersten  (1969) ,  we  define  the  coupling  constant  of  node 
m  by: 


where  ^  ^  is  the  eigenvector  corresponding  to  ^  .  it  is  worth  empha¬ 
sizing  that  the  eigenvalues,  eigenvectors  and  coupling  constants  are 
completely  determined  by  the  material  constants  for  the  thickness  modes  of 
infinite  plates;  these  parameters  are  independent  of  the  boundary  condi¬ 
tions.  The  resonance  frequencies,  in  contrast,  depend  critically  on  the 
boundary  conditions.  The  conditions  that  we  impose  are  that  the  surfaces 
are  traction-free  and  driven  by  a  sinusoidal  voltage: 


j2  *  T2j  -  0  at 

y  =  ±  h 

(20) 

*  *±  *0eiwt 

at  y  =  ±  h 

(21) 

where  we  have  assumed  that  the  plate  is  of  thickness  2h.  We  apply  these 
boundary  conditions  to  our  general  solution  which  is  a  linear  combination 
of  the  modal  eigenvectors.  The  result  is  a  transcendental  equation  whose 
solutions  give  the  allowable  frequencies  of  vibration  for  the  plate.  This 
equation  can  be  expressed  in  the  form  (Tiersten  11969]): 


V'  k  2 
2-r  (n) 


tan  (wh/v(n)) 
w^v(n) 


(22) 
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The  solutions  to  this  equation  yield  the  resonance  frequencies.  This 
equation  cannot  be  solved  in  general  and  demonstrates  the  coupling  of  the 
resonance  frequencies  by  the  boundary  conditions.  In  the  general  case,  the 
frequency  solutions  to  (22)  cannot  be  identified  with  a  certain  mode. 
Fortunately,  we  can  approximate  (22)  by  a  simpler  equation  when  the 
coupling  constants  are  small.  This  approximation  is  very  reasonable  for 
quartz,  where  k^  is  usually  less  than  0.10  (10%) .  If  the  eigenvalues  are 
well  separated  in  addition  to  the  coupling  constants  being  small,  then  the 
resonance  frequencies  can  be  approximated  by  the  solutions  to  the  uncoupled 
equation 


kn2  tan  (wk/vn)  =  l 

(wk/vn 


(23) 


(Ballato  &  Iafrate  [1976]).  Note  that  in  this  approximation  the  modal 
resonances  are  independent.  Hie  solution  to  (23)  can  in  turn  be 
approximated  for  small  kn  by  using  the  Laurent  series  expansion  of  tan  x 
about  x  *  n/2: 


W«m 

nm 


(2m  +  1) 


n  J2 
2  n 


n2  (2m  +  l)2 


(24) 


where  m  =  0,1,...  .  Note  that  the  first  term  is  just  the  solution  to 

wh 

tan  -  =  <v  ,  (25) 

vn 

Hie  correction  term  is  clearly  negligible  for  high  overtones,  i.e.,  for 
large  m  (Bechmann  et  al  [1962]).  We  shall  neglect  this  correction  in  the 
following  and  approximate  the  resonance  frequencies  by  the  solution  to 
(25): 


w. 


nm 


=  (2m  +  1)  —  —  . 


(26) 


Converting  these  angular  frequencies  to  cycles  per  unit  time,  the  resonance 
frequencies  are  given  by 


f 


nm 


n 

(2m  +  1)  — 
4h 


(27) 


13 


The  theory  we  have  sunnarized  so  far  applies  to  a  general,  anisotropic 
piezoelectric  substance;  there  are  no  constraints  placed  on  the  material 
constants.  Hie  most  general  piezoelectric  substance  is  an  arbitrarily 
anisotropic  { triclinic)  crystal  without  center  of  symmetry.  In  this  case, 
there  are  45  independent  material  constants  required  to  specify  the 
elastic,  piezoelectric  and  dielectric  tensors.  However,  we  are  interested 
in  the  properties  of  doubly  rotated  quartz,  so  we  need  not  consider  such  an 
arbitrary  substance.  Quartz  is  a  member  of  the  trigonal,  trapezohedral 
crystal  class  (32,  D^)  (see,  e.g.,  Mason  [1950]).  In  the  crystalline 
coordinate  system,  there  are  only  10  independent  material  constants  (Mason 
[1950],  Tiersten  [1969],  Bechmann  [1956]):  6  elastic  constants,  2 
piezoelectric  constants  and  2  dielectric  constants.  These  constants  have 
been  measured  by  Bechmann  (1956) .  The  elastic,  piezoelectric  and 
dielectric  tensors  for  an  arbitrary  cut  of  quartz  may  be  determined  by 
applying  the  tensor  transformation  rules  to  the  tensors  in  the  crystalline 
coordinate  system.  The  most  general  transformation  between  two  coordinate 
systems  is  given  by  the  Euler  angle  formulation.  (Goldstein  [1950] ) .  Such 
a  transformation  is  specified  by  the  Euler  angles  (  x  ,  M  ,  9  )  (using  the 
convention  of  Goldstein).  We  are,  however,  interested  in  doubly  rotated 
quartz  crystals;  the  plate  orientations  are  specified  by  the  two  angles 
<f>  and  0  relative  to  the  crystal  axes  (see  Figure  2) .  The  case  of  both 
<J>  and  0  being  zero  results  in  the  so-called  Y-plate  (i.e.,  y-axis  normal 
to  plate)  with  material  constants  given  by  those  in  the  crystal  system. 
For  non-zero  <t>  and  0  ,  the  material  constants  are  determined  by  applying 
the  transformation  x  to  the  unrotated  tensors,  where 


x  = 


-sin<*>  cosO 
sin*  sin0 
cos  <t> 


cos*  cose  sin0 

-cos*  sin0  cose 

sin*  ° 


(28) 


If  we  denote  the  tensors  relative  to  the  crystal  axes  by  primes,  the 
elastic,  piezoelectric  and  dielectric  tensors  for  the  doubly  rotated  cuts 
are  given  by: 
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ijkl  =  ximxjnxkoxlpc'nnop 


eijk  =  ^^n^kn^lmn 


(29) 


^ij  xikxjl  £'kl 


Thus,  given  <i>  and  0,  we  calculate  the  material  constants  from  (29).  The 
eigenvalues,  eigenvectors,  coupling  constants  and  approximate  resonance 
frequencies  are  then  calculated  using  these  rotated  tensors. 


The  resonance  frequencies  of  piezoelectric  plates  are  functions  of 
temperature.  The  standard  approach  to  this  temperature  behavior  is 
phenomenological  and  assumes  that  the  plate  is  in  thermal  equilibrium.  We 
then  assume  that  the  equilibrium  resonance  frequency  is  an  analytic 
function  of  temperature  and  expand  it  in  a  Taylor  series  about  a  reference 
state  with  frequency  fQ  and  temperature  TQ: 


OO 


£  ‘co  + 


£(-, 

n=l  \  n 


(n). 


dT 


(n) 


T=T 


) 


(T-T0) 


(30) 


We  may  rewrite  this  in  terms  of  a  fractional  change  in  frequency: 

(T-TQ)n  (31) 


f-fo  Afo  V  ! 

•  i  i  a(n)  f 

- 

Hi 

o 

• 

Hi 

o 

n 

_nj  f  d  T(n) 

T=T 

oJ 

The  standard  notation  for  the  temperature  expansion  coefficients  is  Tf^n^, 
or  explicitly: 


*n  (n)  _  _____ 
f  n  j  tr 


(n)f 


d  T 


(n) 


(32) 


T=T 


Indeed,  we  may  write  the  temperature  expansion  of  any  quantity  g  about  a 


reference  state  (g0,T0)  as: 


a  -  E  T w 


n=l 


g  (T-T0)n  • 


(33) 


It  has  been  empirically  established  that  a  cubic  expansion  suffices  to 


adequately  describe  the  temperature  effects  on  the  quartz  resonances 
(Bechmann  [1956] .  Ballato  [1977] ) .  We  shall  refer  to  the  cubic  frequency 
expansion  as  the  relative  frequency  difference  function,  ^  f/f : 


/\  £  ^ 

-  -  H  Tf  (n)  (T-T0)n  (34) 

f  n=l 


Furthermore,  we  will  assume  that  a  cubic  expansion  suffices  to  describe  the 
temperature  behavior  of  other  material  parameters,  such  as  the  elastic 
constants.  This  is  a  reasonable  approximation  provided  the  temperature 
range  is  small. 


We  shall  now  proceed  to  relate  the  temperature  coefficients  of  the 
relative  frequency  difference  function  to  the  temperature  coefficients  of 
the  material  parameters  for  quartz.  Wfe  shall  do  so  only  for  the  funda¬ 
mental  frequencies.  We  have,  from  (27).  the  modal  frequencies  for  an 
infinite  plate  of  thickness  2k: 


f 


n 


(35) 


We  shall  approximate  the  infinite  plate  by  a  finite  one  of  dimensions  d  and 
1.  This  is  reasonable  provided  d  »  2h  and  1  »  2h.  The  density  is  then 
given  by 


M 

p  s  - 

2hld 


The  modal  frequencies  are  thus: 


Id 


fn  = 


2)M 


(36) 


(37) 


where  M  is  the  plate  mass.  Straightforward  differentiation  gives: 


Note  that  the  temperature  coefficients  and  are  just  the 

coefficients  of  thermal  expansion  in  the  x.  y  and  z  directions.  These  are 


18 


measurable  quantities.  The  temperature  coefficient  of  the  eigenvalue, 
T^n^,  can  be  expressed  in  terms  of  the  material  constants.  The 
characteristic  equation  {16)  yields  a  cubic  for  xn,  which  we  write  as: 

i3-A  i2  +  Bl  +  D=0  (39) 

After  differentiation  and  algebraic  manipulations,  we  obtain: 


T  U>  -11 

\\l 

2  A  T.ID 
n  A 

-  Xn 

B  Tb(1)  -  D  Td(1) 

/  *n\ 

Tx  lx 

n  \  1 

3  xn2- 

2  A  x 

n  +  B 

\ 

Hie  temperature  coefficients  of  A,  B  and  D  are  determined  by  the  eigenvalue 
equation.  We  have: 

R  "  Pu 

+ 

P22  + 

p3 

(41) 

B-  Pu 

P22 

+  \~22 

P33 

+  Pu  P33 

~  P23 

2  _n  2 

I  13 

-  Pu  2 

(42) 

R3  2 

♦P22 

Pl32 

*  P 33  Pl2 

-Pn 

P22  i~33 

-  2 

Pi*  P23  P33 

"2  I 

1  n 

12  1  23 

P* 

(43) 

where 


*2jk2  “  c2jk2  + 


e22je2k2 

^22 


and  for  ease  of  notation,  we  have  written  | jk  instead  of  1 2 jk2 


P 

I2ik2* 


In  a  straightforward  but  tedious  fashion,  the  temperature  coefficients 
of  A.  B  and  C  may  be  found  to  be: 
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T  (1)  = 
AB 


fu  fi 


(1)  +  T  (1)  +  H  P  rp  (1)  + 

+  Tn  +  '  22  >  33  Tn  + 


Pu  t 


(1)  +  T  (1)  .  2  P  2  T  (1) 

p  V  Z  1 23  r  n 

1 33  111  1 23 


.  2  P  2  T  (1)  -  2  P  2  T 

^  *  i3  tp  2  1 12  Tn 


rD(1)  *  D  |  111  T(i^(1)  l~23  i~22  [33  +  2  ^3  T|T^(1)ill  l~23  "R.2  1 13 

+  r 22  Tp  (1)  (ll  2  '^33  Ql  +  2  Ql  Tp  U)  ^2  I  31  ~^23  ^21 
I  22  I  31 

+ 1~33  Tp (i)  ru  -fii  P22 +  2  n.2  Tp a>  Q3  n.2  ■  Qi  Q2) 


The  last  step  in  this  derivation  of  the  first  order  temperature  coefficient 
of  frequency  is  to  relate  the  temperature  coefficients  of  ^  to  those  of 
the  elastic,  piezoelectric  and  dielectric  constants.  This  is  done  using 
equation  ( 44) .  We  shall  assume  that  the  temperature  variation  of  the 
piezoelectric  and  dielectric  constants  can  be  neglected  compared  to  that  of 
the  elastic  constants.  Thus,  from  (44),  we  have  in  this  approximation: 

rjk 

T  (1)  =  -  Tc  (1)  (48) 

|  jk  c2jk2  2jk2 


if 


Let  us  sumrarize  these  results.  We  have  seen  hew  the  first  order  tempera¬ 
ture  coefficient  of  the  relative  frequency  difference  function  may  be 
calculated  from  the  material  constants  and  from  the  first  order  temperature 
coefficients  of  the  elastic  constants  and  the  coefficients  of  thermal 
expansion.  The  steps  in  this  calculation  are  straightforward  but  algebra¬ 
ically  tedious.  The  calculation  of  the  second  and  third  order  temperature 
coefficients  of  frequency  in  terms  of  the  coefficients  of  theriral  expansion 
and  the  elastic  constants  follows  the  same  procedure  and  will  not  be  done 
here.  (Kahan  [1982] ,  Ballato  [1977] ) . 


The  temperature  coefficients  of  the  elastic  constants  and  the  coeffi¬ 
cients  of  thermal  expansion  depend  on  the  orientation  angles  of  the  quartz 
plate  with  respect  to  the  crystalline  axes.  This  clearly  implies  that  the 
frequency  temperature  coefficients  depend  on  these  orientation  angles. 
Thus,  we  need  to  know  hew  to  calculate  the  temperature  coefficients  of  a 
particular  plate  given  those  in  the  crystalline  axes.  Therefore,  we  need 
to  determine  the  transformation  properties  of  the  temperature  coefficients. 
We  will  consider  a  general  tensor  of  order  n,  denoted  by  g-  <  (T) 

1  f  •  •  • 

where  we  have  explicitly  included  the  temperature  dependence.  Now,  for  any 
temperature  T,  the  transformation  between  any  2  coordinate  systems  is 
specified  by: 


3-s  a  (T)  =  xi  k  ...  X.:  ,k  g'k  k  ( 
3l"*,3n  D1'K1  ]n  n  V***  n 


where  x—  is  the  coordinate  transformation  from  the  primed  to  the  unprimed 
axes,  we  therefore  have: 


Now,  the  temperature  expansions  in  the  2  coordinate  systems  may  be  written: 


A  94  i  (T)  =  g;  i  (T0) 

V***  xn  h'*--  xn  ° 


Z  T  (j)  .  (TVT0)i 

j  —1  i  f  .  •  •  i 

1  n 


(51) 


and  for  the  primed  axes: 


A  g'k  k  (T)  =g'k  k  (T  ) 
in  in 


Z  T  ,  (T-T0)j 

j=l  ^  k 

1  n 


(52) 


(In  (51)  and  (52),  there  is  no  implied  suimation.) 

Combining  (50),  (51)  and  (52),  and  equating  like  powers  of  (T-JT0),  we 
obtain  the  following  relation  between  the  temper ature  coefficients  in  the 
two  coordinate  systems: 


T  (j) 

gi  < 

X1 » . . .  in 


1 


^ij/ * • • 


xi  'k 
n  n 


(53) 


Note  that  in  this  formula  there  is  no  implied  summation  over  the  indices 
il,...,in  on  the  right-hand  side.  Using  the  transformation  rule  (53),  the 
temperature  coefficients  of  the  elastic  constants  (a  fourth  rank  tensor) 
for  a  doubly  rotated  plate  can  be  calculated  from  those  in  the  crystalline 
coordinate  system  using  the  coordinate  transformation  (28) .  Similarly,  the 
coefficient  of  thermal  expansion,  a  vector  (i.e.,  first  rank  tensor)  can 
also  be  calculated  from  (53).  The  temperature  coefficients  of  frequency 
can  thus  be  calculated  for  any  doubly  rotated  quartz  plate  from  those  of 
thermal  expansion  and  the  elastic  coefficients  in  the  crystalline  system. 
Indeed,  the  coefficients  for  the  elastic  constants  in  the  crystalline 
system  were  determined  by  Bechmann  (1956)  by  transforming  the  measured 
coefficients  of  various  standard  cuts  of  quartz. 
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4. 


HJMERICftL  RESULTS 


We  shall  new  use  the  theory  of  the  previous  two  sections  to  perforin 
calculations  for  doubly  rotated  quartz  plates.  All  calculations  were 
performed  on  the  Air  Force  Geophysics  Lab's  (APGL)  CDC  6600  computer .  The 
material  constants  and  their  temperature  coefficients  are  those  given  by 
Bechmann  et  al  (1962) . 

We  begin  by  depicting  in  Figure  3  the  values  of  the  coupling  constants 
for  the  three  modes  of  quartz.  The  coupling  constants  were  calculated  from 
(19)  for  the  doubly  rotated  cuts.  Ballato  (1977)  presented  similar  graphs 
using  a  different  definition  of  the  angles  defining  the  doubly  rotated 
cuts.  The  curves  of  Figure  3  are  a  representative  sampling  of  the  coupling 
constant  for  doubly  rotated  cuts.  (Due  to  the  symmetry  of  the  quartz 
crystal  lattice,  we  need  oily  consider  <t>  in  the  range  (0°,60°)  and  0  in 
the  range  (0°,  90°) ) .  We  shall  have  frequent  recourse  to  these  curves  in 
the  following  search  for  cuts  with  desirable  temperature  characteristics; 
for  a  cut  to  be  useful,  it  most  teve  a  reasonable  amount  of  piezoelectric 
coupling. 

The  widespread  applications  c£  quartz  plates  owe  much  to  the  early 
discovery  of  the  AT  cut  ( <t>  *  0°,  0  =  35.0931°) .  From  Figure  3,  we  see  that 
the  AT  cut  has  a  relatively  high  (  ^  9%)  coupling  constant  for  the  C  mode 
and  zero  coupling  for  modes  a  and  b.  In  Figure  4,  we  plot  the  relative 
frequency  difference  function  (mode  C)  about  the  reference  temperature  TQ  = 
25°C  for  the  temperature  range  (-40,80).  The  graph  exhibits  the 
temperature  stability  that  makes  the  AT  cut  the  preferred  choice  in  many 
applications.  This  stability  derives  from  the  flatness  of  the  curve  in  the 
vicinity  of  25°C.  This  flatness  is,  of  course,  a  consequence  of  the 
derivative  vanishing  at  25°,  or  equivalently,  the  first  order  temperature 
coefficient  of  frequency  (for  mode  C)  having  a  zero  at  <t>  ~  0°  and  0  = 
35.0931. 
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Another  of  the  oonmonly  used  cuts  is  the  EfT  cut,  defined  by  <t>  =  60°  and  e  = 
48.9942°,  The  relative  frequency  difference  for  mode  b  is  displayed  in 
Figure  5.  Like  the  AT  cut,  the  BT  cut  owes  its  temperature  stability  to 
the  vanishing  of  the  first  order  temperature  coefficient  of  frequency.  The 
general  shapes  of  the  AT  and  BT  frequency  curves  are,  however,  very 
different.  The  BT  curve  is  parabolic;  it  is  the  second  order  temperature 
coefficient  which  dominates  the  cubic  frequency  expansion.  In  contrast, 
the  AT  curve  is  clearly  of  cubic  form.  These  2  general  forms,  i.e., 
parabolic  and  cubic,  are  representative  of  the  shapes  of  the  relative 
frequency  difference  curves  for  arbitrary  <t>  and  e  . 

The  conmon  feature  of  the  AT  and  BT  cuts  that  rrakes  these  cuts  useful 
is  the  vanishing  of  the  first  order  temperature  coefficients.  Bechmann  et 
al  (1962)  mapped  out  the  locations  of  the  zeroes  in  the  first  order 
temperature  coefficient.  In  Figure  6,  the  loci  of  zeroes  is  displayed  for 
modes  b  and  c  (mode  a  has  no  zeroes) .  This  previous  work  represents  a 
search  for  quartz  cuts  with  stable  temperature  behavior  using  a  very 
"local"  measure  of  sensitivity.  If  the  first  order  temperature  coefficient 
vanishes,  the  changes  in  frequency  due  to  a  small  deviation  in  temperature 
from  25°C  will  also  be  small  (order  of  (T-T)^).  However,  large  temperature 
deviations  may  produce  large  changes  in  frequency,  as  is  evident  for  the  AT 
and  BT  cuts  from  Figures  4  and  5. 

In  many  applications,  the  range  of  temperature  variation  is  too  large 
to  rely  on  a  local  measure  of  sensitivity.  For  these  cases,  a  more 
"global"  measure  is  required,  that  is,  a  measure  which  takes  into  account 
the  variations  over  the  entire  temperature  interval  of  interest.  Further¬ 
more,  modern  electronic  devices  are  capable  of  compensating  for  or 
correcting  small  frequency  changes;  this  is  accomplished  by  following  a 
known  relative  frequency  difference  curve. 

The  first  "global"  measure  that  we  will  consider  uses  an  integral 
measure  of  temperature  sensitivity.  Specifically,  the  measure  is  the 
integral  of  the  absolute  value  of  the  relative  frequency  difference  from 
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AN6LE  OF  ROTATION  -  THETA 


r 


-40°  to  80°,  i.e. , 

(  <t> ,  e  ) 


dT 


(54) 


where  the  notation  explicity  displays  the  dependence  of  the  measure  on  the 
cut  angles  <t>  and  6  .  In  Figure  7,  we  present  the  results  of  the  calcula¬ 
tions  of  (0,0)  for  the  three  modes  of  quartz  plates.  These  curves 
give  us  a  rough  indication  of  the  location  and  magnitudes  of  the  minima  in 
temperature  sensitivity.  The  minina  for  mode  a  are  seen  to  be  very  weak 
and  will  not  be  considered  further.  The  minima  for  modes  b  and  c  are  of 
reasonable  order  (i.e.,  relative  magnitude  of  approximately  100)  and 

warrant  further  investigation.  The  actual  locations  and  magnitudes  of  the 
b  and  c  mode  minima  were  calculated  by  minimizing  (0,0)  with  respect 
to  0  for  fixed  <t> .  The  results,  for  <t>  in  increments  of  1°,  are  given  in 
Figure  8  for  mode  c  and  Figure  9  for  mode  b.  Hie  absolute  minima  for 
mode  c  occurs  for  <t>=  0Q,  while  for  mode  b,  it  occurs  for  0  =  60°.  It  is  no 
coincidence  that  these  are  the  0 -values  of  the  AT  and  BT  cuts, 

respectively.  However,  the  0  values  of  the  minima  do  not  agree  exactly 
with  those  of  the  AT  and  BT  cuts;  the  minima  do  not  have  vanishing  first 
order  temperature  coefficients.  For  mode  c,  the  minima  occurs  at  0  =  0°, 
0  =  35.13C43;  for  mode  b.  it  occurs  at  0=  60°,  0  =  49.3833°.  The 

relative  frequency  difference  curves  of  these  minima  are  given  in  Figures  4 
and  5  alono  with  the  corresponding  standard  cut.  The  minima  for  modes  b 
and  c  thus  occur  in  very  close  proximity  to  the  zeroes  of  the  first  order 
temperature  coefficients.  This  is,  in  fact,  true  of  all  the  minima 
displayed  in  Figures  8  and  9.  Indeed,  if  we  were  to  plot  the  0  and  0 

angles  of  the  minima  in  Figures  8  and  9,  we  would  produce  a  curve  very 
similar  to  Figure  6,  the  loci  of  the  zeroes  of  the  first  order  temperature 
coefficient. 


We  have  seen  that  the  least  sensitive  cut,  with  respect  to  an  integral 
measure  is  given  by  mode  c.  0  =  0°,  0  =  35.130^.  The  c  mode  minima  is 
superior  by  more  than  a  factor  of  10  to  the  b  mode  minima.  In  fact,  from 
Figure  8,  we  see  that  any  c  mode  minima  with  0<50°  is  superior  to  the  b 
mode  minima.  However,  we  must  use  more  than  the  magnitude  of  the  minimum 
to  compere  different  cuts  effectively.  The  goal  here  is  to  find  cuts  that 
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Figure  7:  INTEGRAL  MEASURE  OF  TEMPERATURE  SENSITIVITY  (continued) 
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Figure  7:  INTEGRAL  MEASURE  OF  TEMPERATURE  SENSITIVITY  (continued) 
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Figure  7:  INTEGRAL  MEASURE  OF  TEMPERATURE  SENSITIVITY  (continued) 
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Figure  7:  INTEGRAL  MEASURE  OF  TEMPERATURE  SENSITIVITY  (continued) 
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Figure  7:  INTEGRAL  MEASURE  OF  TEMPERATURE  SENSITIVITY  (continued) 


Figure  9:  INTEGRAL  MINIMA  FOR  C  MODE 
(MINIMA  WITH  RESPECT  TO  9  )  . 
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are  useful  and  so  various  practical  criteria  must  also  be  considered.  The 
first  such  criterion  we  will  treat  is  angle  stability:  for  a  minimum  to  be 
useful,  the  temperature  sensitivity  must  not  change  greatly  if  the  angle  is 
changed  slightly.  Such  angular  deviations  or  errors  are  likely  to  occur  in 
the  manufacturing  process.  To  evaluate  this  angular  stability  numerically, 
we  will  consider  the  effect  of  errors  in  the  0  value  of  10'  (0.167°).  For 
the  minima  in  Figures  8  and  9,  the  integral  measure  was  also  calculated  for 
the  same  <*>  value ,  but  with  0  changed  by  +10'  and  -10'.  The  upper  curves 
in  Figures  8  and  9  represent  the  larger  of  the  integrals  for  the  +10'  and 
-10'  cases.  It  is  thus  a  reasonable  measure  of  angular  stability.  The 
deficiencies  of  the  <t>  =  0°,  mode  c  minimum  are  apparent:  it  is  highly 
sensitive  to  stall  changes  in  angle.  Indeed,  a  10'  change  in  0  produces  a 
relative  change  in  the  integral  measure  of  several  hundred  percent.  This 
effect  of  snail  0  errors  on  the  relative  frequency  difference  is  demon¬ 
strated  in  Figure  10.  This  angular  instability  is  also  demonstrated  by  the 
very  sharp  minima  for  mode  c  at  0=0°  given  in  Figure  7.  In  contrast,  the 
minima  for  <t>  =  60°  and  mode  b  is  seen  in  Figure  7  to  be  somewhat  rounded 
and  not  as  abrupt.  This  is  borne  out  by  Figure  9,  where  a  10'  deviation  in 
0  produces  a  relative  change  of  about  2%  in  the  temperature  measure. 
Thus,  the  b  mode  minimum  is  very  stable.  This  stability  is  exhibited  by 
the  relative  frequency  difference  curves  in  Figure  11.  Because  of  this 
extreme  angular  instability,  the  <t>  -  0°,  c  mode  minimum  cannot  be 
practically  considered.  However,  as  mentioned  previously,  there  are  other 
c  mode  minima  which  have  analler  minima  than  the  b  mode.  Taking  into 
account  angular  stability,  inspection  of  Figure  8  reveals  two  candidates 
for  useful  c  mode  cuts:  =  27°  (  6  =  34.064^)  and  <t>  =  42°  (  0  = 

33.2906°) .  The  effect  of  theta  errors  on  the  relative  frequency  difference 
for  these  cuts  is  depicted  in  Figures  12  and  13. 

A  second  practical  criterion  is,  as  mentioned  earlier,  the 
piezoelectric  coupling  constant.  The  mode  b  minimum  (<P=  60°,  0  = 

49.3833°)  is  very  close  to  the  standard  BT  cut  (  <fi  =  60°,  0  =  48.9942°). 

Thus,  we  would  expect  this  minimum  to  satisfy  most  practical  criteria  with 
regard  to  the  coupling  constant.  For  this  reason,  we  will  restrict  the 
calculations  of  coupling  to  the  mode  c  minima,  especially  the  two  mentioned 
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above.  In  Figure  14,  we  display  the  results  of  the  calculation  of  the 
coupling  constant  for  mode  c  along  the  mode  c  minimum  curve  of  Figure  8. 
Hus  graph  shows  that  for  <t>  =  27°,  the  minimum  has  a  coupling  constant  of 
about  3.7%,  while  it  is  approximately  2.2%  for  the  <P  =  42°  minimum.  By 
comparison#  the  mode  b  minimum  (  <t>=  60°)  has  larger  piezoelectric  coupling, 
approximately  5.5%  (see  Figure  3). 

Hie  last  practical  criterion  we  will  consider  pertains  to  the  separa¬ 
tion  in  frequency  of  the  three  modes.  From  equation  (27)  we  have  that  the 
fundamental  frequencies  of  the  three  modes  may  be  written  as: 

1 

f"=  w? 

In  this  form,  it  is  clear  that  we  may  treat  the  separation  of  eigenvalues 
instead  of  frequency.  In  Figure  15,  the  eigenvalues  of  all  3  inodes  are 
plotted  along  the  mode  c  minimum  curve  of  Figure  8.  For  both  of  our 
prospective  c-mode  minima,  the  c-mode  eigenvalue  is  well  separated  from 
those  of  the  other  two  modes.  This  criterion  is  not  really  applicable  to 
our  b-mode  minimum  because,  for  <t>  =  60°  and  0  -  49.3833°,  the  piezoelectric 
coupling  vanishes  for  inodes  a  and  c  (see  Figure  3) .  There  is  thus  no 
difficulty  in  this  case  due  to  the  interference  of  the  plate  vibrations  of 
the  other  inodes.  This  is  an  advantage  of  the  mode  b  minima  over  our  2  mode 
c  minima,  which  have  piezoelectric  vibrations  for  all  three  modes  present. 

Let  us  summarize  our  results  for  the  integral  measure  of  temperature 
sensitivity.  We  have  located  four  cuts  which  minimize  this  sensitivity  and 
have  outlined  sane  of  the  advantages  and  disadvantages  of  these  cuts.  The 
absolute  minimum  (0=0°,  mode  c)  is  extremely  sensitive  to  errors  in  the 
cut  angles  and  this  instability  renders  precision  applications  very 
difficult.  Of  the  remaining  three  cuts,  the  weakest  minimum  (  <j>  =  60°, 
mode  b)  possesses  the  greatest  angular  stability,  the  largest  coupling 
constant,  and  has  only  one  mode  present.  However,  the  range  of  variation 
of  the  relative  frequency  difference  curve  may  be  too  large  (i.e.,  the 
minimum  too  weak)  to  be  practical. 
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Next,  we  shall  consider  another  measure  of  sensitivity.  Define  the 
range  measure  by: 


M2(4>,0) 


A* 


nax  — 
f 


(56) 


where  the  maximum  and  minimum  are  over  the  temperature  range  (-40°,  80°) . 
As  for  the  integral  measure,  we  calculate  ^  (<£  ,e  )  for  each  mode  for  <t>  in 
increments  of  10°  and  0  in  the  range  0°  to  90°.  The  results  are  displayed 
in  Figure  16.  It  is  instructive  to  compare  Figures  16  and  7.  It  is  seen 
that  the  shapes  of  the  curves  are  virtually  identical.  Similar  calcula¬ 
tions  were  performed  for  a  third  measure  of  temperature  sensitivity  using  a 
slope 


M3(  <b  ,  0  )  =  max 


d  A f 

dT  f 


(57) 


The  curves  for  this  measure  will  be  emitted  since  they  are  indistinguish¬ 
able  from  those  of  the  other  two  measures.  For  all  three  measures,  the 
minima  for  mode  a  are  weak  and  are  not  treated  further.  Using  the  measures 
M2  and  M3,  curves  of  the  minima  for  fixed  0  were  calculated  for  modes  b 
and  c  for  4>  in  increments  of  1°.  The  resulting  curves  are  not  shewn  since 
they  ar^  identical  in  shape  to  Figures  8  and  9.  It  follows  that  the 
overall  features  of  the  temperature  sensitivity  using  Mj  and  M3  are  the 
same  as  for  M^.  in  all  three  measures,  we  can  identify  three  cuts  (  <j>  = 
0°,  27°  and  42°)  which  minimize  the  mode  c  sensitivity  and  one  (  <f>=  60°) 
which  minimizes  that  of  mode  b.  We  summarize  the  magnitudes  and  locations 
of  these  minima  in  Table  1.  It  can  be  seen  that  the  locations  of  the 
minima  for  different  measures  differ  only  very  slightly.  Indeed,  the 
angular  locations  of  these  minima  deviate  only  slightly  from  the  locations 
of  zero  first  order  temperature  coefficients.  It  can  also  be  shown  that 
the  features  of  these  minima  with  regard  to  angular  stability,  coupling  and 
frequency  separation  are  the  same  for  all  three  measures.  The  relative 
advantages  and  disadvantages  of  these  4  minima  are  the  same  for  all  three 
measures  and  are  as  stated  earlier. 
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Figure  16:  RANGE  MEASURE  OF  TEMPERATURE  SENSITIVITY 
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Fic^re  16:  RANGE  MEASURE  OF  TEMPERATURE  SENSITIVITY  (continued) 
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Figure  16:  RANGE  MEASURE  OF  TEMPERATURE  SENSITIVITY,  (continued) 
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Fic^re  16:  RANGE  MEASURE  OF  TEMPERATURE  SENSITIVITY  (continued) 
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Figure  16:  RANGE  MEASURE  OF  TEMPERATURE  SENSITIVITY  (continued) 
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Figure  16:  RANGE  MEASURE  OF  TEMPERATURE  SENSITIVITY  (continued) 
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Figure  16:  RANGE  fEASURE  OF  TEMPERATURE  SENSITIVITY  (continued) 
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Figure  16:  RANGE  MEASURE  OF  TEMPERATURE  SENSITIVITY  (continued) 
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Figure  16:  RANGE  f-EASURE  OF  TEMPERATURE  SENSITIVITY  (continued) 


PHI  =  40  MODE  =  3 


I  OG  (RANGE) 


PHI  «  50  MODE  *  3 


Figure  16:  RANGE  MEASURE  OF  TEMPERATURE  SENSITIVITY  (continued) 
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Figure  16:  RANGE  MEASURE  OF  TEHFERATORE  SENSITIVITY  (continued) 


TABLE  1:  Minina  Using  Measures  M^,  M2  and  M3 


MODE 

1 

e 

1 

M^  min 

0  Mx 

M2  min 

e  m2 

m3 

© 

C 

0 

35.1304  . 4035xl0-3 

35.1720  ,1181xl0-4 

35.3242 

C 

27 

34.0644  .2246xl0-2 

34.0216  .5897X10"4 

34.1258 

c 

42 

33.2906  .1564xl0~2 

33.1912  . 4102x1 0“4 

33.1697 

b 

60 

49.3832  .5385xl0-2 

48.9622  .1393xl0~3 

48.5290 

.3782X10-6 
.196 4x1 0-5 
.136 9x1 0-5 
. 4495xl0-5 
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The  three  measures  that  we  have  investigated  all  measure  the  sensi¬ 
tivity  of  a  single  mode.  A  slight  generalization  of  these  methods  uses  two 
modes,  that  is.  measures  the  temperature  sensitivity  of  the  difference  in 
frequency  between  two  modes.  Let  mode  i  have  temperature  coefficient  (j 
=  1,2,3)  and  mode  j  have  coefficient  b..  (j  =  1,2,3).  By  definition,  we 
have: 


J 

fi(T)  =  fi(T0)  j\  +  Lak  (T  -  T0)k  j 
fj(T)  -fj(T0)£l  +  Z  bk  (T-T0)kj 


Define  the  frequency  difference  between  modes  i  and  j  by: 


Sij(T)  =  f i(T)  -  f j  (T) 


Substitution  of  (58)  and  (59)  gives: 


sij(T)  =  Si;j(T0)  +  £  akfi(T0)  -bkfj(T0)  (T  ~  TQ) 
k=l  '  ’ 


Hie  fractional  change  in  the  frequency  difference  is  thus: 


ASy  Sij{T)  -S^IV 


S«(T„) 


where  we  have  defined  the  weights: 


fi(Tp) 

fi(T0)  “  fj(T0) 


(64) 


Wj  = 


£j  ffp) 

fi<T0>  "  fj<To) 


(65) 


These  weights  are  dimensionless  quantities  which  can  be  written  in  terms  of 
the  eigenvalues.  For  the  fundamental  frequencies  of  modes  i  and  j,  substi¬ 
tution  of  (55)  into  (64)  and  (65)  gives: 


Wi  = 


(66) 


Wj  = 


(67) 


These  can  also  be  written  in  terms  of  the  modal  velocities  as: 

v. 


Wi  = 


v.  -  y 
1  J 


(68) 


Wj  =  - ^ -  (69) 

Vi  "  vj 

Equation  (63) ,  with  the  weights  defined  by  (66)  and  (67) ,  was  used  to 
investigate  the  temperature  sensitivity.  Note  that  the  quantity  S^j(T)  is 
a  directly  measurable  quantity.  The  sensitivity  was  measured  using  the 
integral  of  the  absolute  value  of  (63)  from  -40°  to  80°.  The  calculations 
and  procedure  were  identical  to  those  for  the  single  mode  measures.  The 
results  are  summarized  graphically  in  Figure  17.  These  graphs  should  be 
compared  to  the  single  mode  data  of  Figure  7.  This  comparison  shows  that 
all  the  minima  of  the  double  mode  measure  are  very  weak.  None  of  the 
minima  of  Figure  17  is  less  than  0.01.  In  contrast,  our  single  mode 
integral  measure  has  minima  much  less  than  this  value  (see  also  Table  1) . 
The  double  mode  measure  will  not  be  investigated  further  since  it  is 
inferior  to  the  single  mode  measures. 
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Fi^ire  17:  DIFFERENCE  JCASURE  OF  TEMPERATURE  SENSITIVITY 
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Figure  17: 
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Figure  17;  DIFFERENCE  MEASURE  OF  TEMPERATURE  SENSITIVITY  (continued) 
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Figure  17:  DIFFERENCE  JBASURE  OF  TEMPERATURE  SENSITIVITY  (continued) 
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Figure  17: 
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Figure  17:  DIFFERENCE  JEASURE  OF  TEMPERATURE  SENSITIVITY  (continued) 
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Previous  work  on  the  temperature  sensitivity  of  quartz  resonances  used 
a  local  measure  of  temper ature  sensitivity,  namely,  the  zeroes  of  the  first 
order  temperature  coefficients.  We  have  presented  results  of  calculations 
using  other  measures  of  sensitivity  in  the  temperature  range  (-40°,  80°) 
and  have  found  that  the  minima  for  fixed  (j>  occur  in  close  proximity  to  the 
locus  of  zeroes  of  the  first  order  temperature  coefficients.  Indeed,  we 
have  found  that  the  locations  of  the  minima  vary  only  slightly  between  the 
different  measures;  the  locations  of  the  minima  for  the  integral  measure, 
the  range  measure  and  the  slope  measure  are  all  in  close  proximity  to  each 
other  and  to  the  zeroes  of  the  first  order  temper  ature  coefficient.  It 
seems  likely  that  this  characteristic  is  attributable  to  the  cubic  nature 
of  the  relative  frequency  difference  expansion.  Nonetheless,  our  results 
do  go  beyond  an  analysis  of  the  zeroes  of  the  first  order  temperature 
coefficient.  Specifically,  the  measures  we  have  used  allow  us  to  compare 
the  temperature  sensitivity  of  the  various  minima  according  to  their 
magnitude.  In  addition,  we  have  used  three  other  practical  criteria  in 
comparing  these  doubly  rotated  cuts;  angular  stability,  separation  of  modal 
frequencies  and  coupling  strengths.  In  this  fashion,  we  have  presented 
four  cuts  with  relatively  stable  temperature  behavior.  None  of  these  four 
cuts  is  clearly  superior  to  the  others,  but  each  has  its  advantages  and 
disadvantages. 
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